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ABSTRACT 



The strength of Titan's methane cycle, as measured by precipitation and 
evaporation, is key to interpreting fluvial erosion and other indicators of the 
surface-atmosphere exchange of liqui ds. But th e mechanisms behind the 



occurrence of large cloud outbursts (jSchaller et al.l 120061 ) and precipitation 
([Turtle et al.ll2009l . l2011bl ) on Titan have been disputed. A gobal- and annual- 
mean estimate of surface fluxes indicated only 1% of the insolation, or ~0.0 4 
W/m 2 , is exchanged as sensible and/or latent fluxes (jMcKav et al.l Il991 



Since these fluxes are responsible for driving atmospheric convection, it has 
been argued that moist convection should be quite rare and precipitation even 
rarer, eve n if evaporation glo bally dominates the surface-atmosphere energy 
exchange ([Griffith et al.ll2008l ). In contrast, climate simulations that allow at- 
mospheric motion indicate a robust methane cyc le with substantial cloud for- 



mation and /or precipitation (IRannou et al 



2008 



Mitchell et al. 



2009, 



2011 



2006; 



Schneider et al 



Mitchell et al. 



2006; 



Mitchell 



20JJ). We argue the top-of- 



atmosphere radiative imbalance - a readily observable quantity - is diagnostic 
of horizontal heat transport by Titan's atmosphere, and thus constrains the 
strength of the methane cycle. Simple calculations show the top-of-atmosphere 
radiative imbalance is ~0.5-l W/m 2 in Titan's equatorial region, which implies 
2-3 MW of latitudinal heat transport by the atmosphere. Our simulation of 
Titan's climate suggests this transport may occur primarily as latent heat, with 
net evaporation at the equator and net accumulation at higher latitudes. Thus 
the methane cycle could be 10-20 times previous estimates. Opposing seasonal 
transport at solstices, compensation by sensible heat transport, and focusing of 
precipitation by large-scale dynamics could further enhance the local, instan- 
taneous strength of Titan's methane cycle by a factor of several. 

Titan's transport-dominated methane cycle is fundamentally distinct from the hy- 
drological cycle of Earth's tropical zone. Earth's hydrology is primarily driven by local 
radiative imbalance at the surface that amounts to ~115 W/m 2 of evaporation occur- 
ring over the tropics. An additional ~40 W/m 2 of evaporation results from horizontal 
heat transport away from the equator by the motion of the atmosphere-ocean system 
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(jTrenberth fc Stepaniakl l2003al) . Following conv entions from the literature on Earth's cli- 
mate energetics (ITrenberth Stepaniakl l2003al ) . we identify the role of energy transport 
by the atmosphere in Titan's global energetics. Vertically integrated quantities are de- 
fined M = J °° pMdz = ± J Ps Mdp. The total energy of the atmosphere, after vertical 
integration, is 



A E = k + c p T + <S> s + Lq 
« D E + L E = M E , 



(1) 
(2) 



with kinetic energy K E = k (assumed to be negligible), internal energy c p T, potential 
energy $ s (which combine to give dry static energy D E = c p T + $ s ) and latent energy 
L E = Lq. The atmospheric energy is very nearly equivalent to the vertical integral of the 
moist static energy, M E = D E + L E = h, with h = c p T + gz + Lq. 



The thermodynamic equation, 

dD E 
dt 

dL E 
dt 



and the moisture equation, 



-v-f De + q 1 



-V • F Le - Q 2 



combine to give the moist static energy equation 

dM E 



dt 



V ■ F Me + Qi 



(3) 



(4) 



(5) 



where 



(6) 
(7) 
(8) 

(9) 
(10) 

(11) 

Q\ is the column-integrated diabatic heating in W/m 2 . Atmospheric heating can arise from 
radiative imbalances at the top-of-atmosphere, Rt, and surface, Rs, from surface sensible 
heat fluxes, H$, or the release of latent heat by precipitation, LP. The column-integrated 
latent heating, Q2, is the difference of precipitation and evaporation, L(P — E). 



Fd e 


= vepT + v$ 


f l e 


= vLq 


Fm e 


= F De + F Le 


Q2 


= L(P-E) 


Qi 


= R T -R S + H S + LP 




= R T + F S + Q 2 . 
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In a steady-state, there is a balance between column heating, Q\ — Q2, and horizontal 
divergence of moist static energy, V • Fm e - In an atmosphere with no horizontal heat 
transport (as is true in single-column, radiative-convective models), this balance reduces 
to Qi = Q2, which further reduces to F s = Rt if the atmosphere does not store heat. 
In this case, the net surface flux F s = LE S + H s — R s is equal and opposite to the top- 
of-atmosphere radiative imbalance Rt = S — OLR, i.e., a heat source/sink at the top of 
the column must be offset locally by a sink/source at the bottom^ Global equilibrium 
requires the average of Rt to be zero, therefore on average the surface radiative imbalance 
R s is offset locally by turbulent exchange of latent heat LE S and sensible heat H s with the 
atmosphere. Based on these arguments, a one-di mensional radiative-convective model of 



Titan with R s ~ 0.04 W/m ([McKay et al.lll99ll ) is often invoked as a constraint on the 



evapo rative surface flux, L E S < R s , a severe limitation on the strength of Titan's methane 



cycle ([Griffith et al.ll2008l ). We argue that horizontal heat transport by the atmosphere 



may dramatically increase this limit. 

In order to derive a constraint on the strength of Titan's methane cycle in the presence 
of meridional heat transport by the atmosphere, first we observe Tita n's OLR is ess entially 



"flat" (independent of latitude) and "steady" (independent of time) ([Li et al.ll201ll . Figure 
[Th,). Time-independence results from the long radiative cooling time of Titan's middle 
troposphere, where the atmosphere radiates to space. However, long thermal times do not 
account for a flat OLR. Indeed without atmospheric heat transport, Titan's OLR would 
exactly balance the annual-mean insolation which is strongly peaked at the equator and 
decreases poleward (Figure [[b, solid line). Instead, Titan's flat OLR (Figures Hk.,b) is 
direct evidence for poleward energy transport by the atmosphere (Figure [2]). Since the 
insolation S is a function of both space and time and we require global energy balance on 
long timescales (i.e., constant Ae), the top-of-atmosphere radiative imbalance Rt ~ S — S a 
with S the global- and time- mean insolation (Figure [TH, dashed line). In this case, Rt is 
offset either by surface energy storage or by horizontal transport of moist static energy. The 
former option is ruled out by the solid, low-heat-capacity surface of Titan which prevents 
energy from being transported or stored for long periods. Thus, we infer a divergent flux 
of moist static energy from regions of positive Rt to regions of negative Rt (Figures [TH 
and [2b). 

The implied heat flux from Rt is in general a combination of latent and dry static 
energy (Figure [2]) , however we argue the imbalance itself provides a rough estimate of the 
magnitude of column cooling by latent heat transport out of the tropics (the "drying rate" 



1 Our definition of S is the net solar flux, and includes the reduction due to local albedo. 
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do to evaporation), 



Q2 



LE ~ 0.5 



W/m 2 



7 cm/yr / 

which is in fact the value obtained in our simulation (Figure [2f). 



(12) 



Our arguments thus far are based on the assumptions of top-of-atmosphere radia- 
tive imbalance, surface energy balance, and the dominance of poleward latent heat trans- 
port. Further insight into Titan's methane cycle can be gained through global circulation 
model (GCM) simulations. In the annual mean, the atmospheric heat transport implied 
by the TOA radiative imbalance is from the equator to higher latitudes. However, simu 



lations flRannou et al.ll2006l ; iMitchell et al. 



Schneider et al 



Turtle et al. 



2006 



Mitchell 



2008; 



20121 ) and cloud observations (jBrown et al.l l2010l ; iRodriguez et al. 



Mitchell et al 



2009, 



2011 



2011 



201 lal ) indicate the transport seasonally reverses from a globally northward 



phase to a southward phase, with large cancellation between the two phases in the annual 
mean. A compensation between dry-static and latent transports may further enhance the 
strength of Titan's methane cycle, as we now demonstrate. 

In the remainder of this paper, we substantiate the above arguments and draw atten- 
tion to differences between Titan's methane cycle and Ear th's hydrological cy cle by com- 
paring GCM simulations of Titan with NCEP reanalysis (jKalnav et al.lll996l ) for Earth. 
We focus our attention on a "moist" version of Titan, i.e., one with an unlimited and pure 
supply of surface methane, in order to establish the potential strength of the methane cycle. 
We then discuss the effects a limited surface supply and the presence of liquid hydrocarbon 
mixtures would have on our results. 

Figure [2] shows a br eakdown of the energetics of our GCM simulation of Titan reported 



m 



Mitchell et al.l (|201ll ). The flux of moist static energy seasonally reverses with peak 



magnitudes of ~11 W/m 2 (|2h.). In the annual mean, there is flux divergence at the equator 
and convergence in mid-latitudes, with peak fluxes of ~2 MW (j2jo). The annual-mean net 
column heating Q\ — Q2 (top-right panel) is positive at the equator and negative at high 
latitudes, indicating an excess of ~0.5 W/m 2 insolation at the equator and an excess of 
> 0.5 W/m 2 OLR near the poles. The divergence of moist static energy offsets the net 
column heating. 

At this point, we encounter the first feature of Titan's climate that distinguishes it 
from Earth's tropical climate. When we break the components of the moist static energy 
transport into latent (Figure [2]i,e,f) and dry static (Figure [2g,h,i) transport, we find that 
the annual- mean heat flux divergence at low latitudes is dominated by latent energy fluxes. 
As a result, the equatorial region experiences net drying (negative Q2] Figure Ef) of ~0.5 
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W/m 2 , while mid-latitudes have net accumulation of a similar magnitude, an effect already 
well documented in several Titan c l imate models (Ranno u et alJl200a : iMitchell et al .1 12006: 

2OI2J). This confirms our earlier 



Mitchell 



2008 



Mitchell et al. 



2009 



2011 



Schneider et al 



assumption, that the magnitude of the equatorial drying rate is set by the horizontal latent 
energy transport required to keep Titan's OLR "flat" . 

A second distinguishing feature of Titan's climate follows by dividing the annual- 
mean latent energy transport into components due to the time-mean and variations as 
Le = Le + L' E , such that for instance Fl<; = Lvq + Lv'q', we find that the majority of 

the energy transport occurs due to the time-dependent seasonal cycle (Lv'q'; dashed line 
in Figure [2b) rather than by the annual-mean component (Lvq; dotted line). In fact, the 
mean component is convergent at the equator, but the seasonal divergence dominates. In 
Earth's tropics, by contrast, the a tmosphere accomplishes moist static energy transport 
primarily by the mean component (jTrenberth &: Stepania k 2003b). 



Our analysis reveals a third important feature of Titan's atmospheric energetics: a 
large compensation between dry-static and latent energy transports during solsticial con- 
ditions (Figure [2b,d,g). In fact, instantaneous latent energy fluxes oppose the moist static 
energy flux (compare Figures [2b, and d), implying that moisture is being fluxed towards 
the summer pole. Instantaneous dry-static transported from the summer to the winter 
hemisphere carries is roughly double the latent transport (Figure [2g), and in the annual 
mean is moderately convergent at the equator (Figure [2h, solid line) . This compensation i s 
also a feature of Earth's annual-mean tropical circulation (jTrenberth &; Stepaniakll2003bl ). 
but in contrast Earth's quasi-steady tropical circulation converges moisture to the equa- 
tor and diverges dry-static energy towards the poles. Based on our simulations, it seems 
plausible that Titan's seasonal cycle enhances the instantaneous methane transport and 
precipitation by at least a factor of five, or ~10 W/m 2 . 

The surface radiative flux imbalance is the energy available to drive turbulent exchange 
of sensible (dry-static) and latent heat flux es, and is often cited as an upper-limit to the 
evaporative flux in Titan's climate system ([Griffith et alJ 120081 ; IWilliams et al.l 12012 ). A 
comparison between annual-mean surface radiative imbalance (Figure [3]) from our "dy- 
namic model" with horizontal heat transport (solid line) and in radiative-convective equi- 
librium with the annual- and global- mean insolation (thin gray line) demonstrates the 
influence of atmospheric circulation on Titan's methane cycle and Earth's hydrological cy- 

For Titan, the surface radiative imbalance in the dynamic model is everywhere larger 



2 In the case of Earth, radiative-convective equilibrium was calculated with a single-column version of 
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than the radiative-convective model surface imbalance. The modeled instantaneous fluxes 
can be a factor of several larger than this (Figure [4]) primarily due to the strong seasonal 
cycle in the TOA radiative imbalance (Figure [1]). By contrast, Earth's surface radiative 
imbalance is dominated by "local" heating that does not involve latitudinal heat transport 
(Figure [3b, thin gray line). Comparing with NCEP fluxes, we see that transport accounts 
for a ~30% enhancement of radiative imbalance at the equator (solid line). This dominance 
of "local" energetics in Earth's climate may be the origin of misconceptions about Titan's 
climate that lead to large underestimates of the potential strength of the methane cycle. 
To reiterate, Titan's "flat" OLR (Figure [TJ is direct evidence for 2-3 MW of latitudinal 
heat transport by the atmosphere (Figure [2b), our numerical experiments suggest much of 
it could be transported as latent energy (Figure [2b), and therefore Titan's methane cycle 
has the potential to be much stronger than is estimated based on global average surface 
energetics. 

There exist a few mechanisms not accounted for in our analysis that could throttle 
Titan's methane cycle. Firstly a nd most important ly, there is apparently a limited supply 
of liquid methane at t he surface (iLqrenz et al.1120081 ). Although precipitation does occur in 
the equatorial region dTurtle et alJboilbh. dry conditions are expected to persist there in 



the present climate (iRannqu et al 



2009 



Schneider et al 



2006 



Mitchell et al.ll2009 : lMitchellll2008l : Mitchell et al. 



20121 ). Sensitivity experiments with our model demonstrate that as 



the global methane supply is reduced, the latent heat transport in the e quatorial regio n 
reduces in-kind and dry static transport increases to offset the reduction (jMitchelll 120081 ) . 
Secondly, liquid ethane may be well-mixed in liquid methane at the surface, and if so it 
would suppress the saturation vapor pressure of methane at a level proportional to the 
unknown mixing ratio of the two hydrocarbons in surface reservoirs. We tested this effect 
in our two-dimen sional climate simulations by artificially suppressing the evaporation of 
surface methane ([Mitchell et al .1 120061 . 120091 ) . Our results indicate that while there is a 
substantial reduction in the annual-mean evaporation rate relative to the "moist" case, the 
instanta neous latent h eat transport at solstices is only reduced by a factor of 2-3, i.e., to 
7 MW (jMitchelll 120081 ) . Local evaporation and accumulation in the annual- mean in this 
scenario are of order 0.5-1 W/m 2 , similar to the "moist" case (Figure Gf). 

In summary, Titan's nearly constant OLR, seasonal lower tropospheric circulation, and 
compensation between dry-static and latent energy transports allow a vigorous methane 
cycle despite weak, global-mean radiative forcing. Titan's storms are triggered in large- 



the CCM3 model forced by the annual- and global-mean insolation. NCEP reanalysis is used in lieu of a 
"dynamic model". 
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scale updrafts (IMitchell et al.l 12 006; 



20121 ) and waves (IMitchell et al 



precipitation events (jTurtle et al 



Mitchell 



2011 



2008; 



Mitchell et al 



2009; 



Schneider et al. 



2009 



Large cloud outbursts (jSchaller et al.ll2006l ) and 



2011bl ) occur as the result of the global transport 



and focusing power of Titan's atmospheric circulation, with significant implications for 
the interpretation of surface features ascribed to the presence of flowing and evaporating 
liquids. 
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Fig. 1. — Titan's top-of-atmosphere (TOA) radiative fluxes are show separately as outgoing 
longwave radiation (OLR) and insolation (S) as a function of seasons in (a & c) on the same 
color scale, and as annual- means or instantaneously as indicated in (b &: c). Titan's OLR 
(a) is nearly independent of latitude and time, while the insolation (c) varies substantially 
with season. In the annual mean, insolation is strongly peaked at and symmetric about 
the equator (solid line). Since the OLR is flat, there is a net imbalance in TOA radiation 
(c), with ~0.3 W/m 2 of heating at the equator and ~0.5 W/m 2 cooling at high latitudes. 
The insolation at southern summer solstice (dotted line) is peaked at the summer pole and 
gives rise to a net TOA imbalance of up to 2 W/m 2 . 
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Fig. 2. — Our model of Titan's climate dynamics as a function of latitude and season (a) 
indicates peak transport of 11 MW. In the annual- mean (b) the transport needed is ~2 
MW to offset the TOA radiative imbalance (c & Figure [1]). Time-dependent (dashed line in 
b) and mean (dotted line) transport both contribute to the moist static energy transport. 
The latent heat transport (d) peaks after solstices at 11 MW, and is offset by dry-static 
energy transport (g) in the opposing direction. The total moist static energy transport 
in the equatorial region (b) is dominated by latent heat transport (e) while dry-static 
transport is nearly negligible. Divergent latent energy transport in the equatorial region 
(e) is balanced by net evaporation near the equator and accumulation in mid-latitudes (f). 
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ta) Titan surface energetics (b) Earth surface energetics 




-0.4 -0.2 0.0 0.2 0.4 -200 -100 100 200 

Flux [W/m 2 ] Flux [W/m 2 ] 



Fig. 3. — Annual-mean net surface radiative fluxes, which are responsible for driving 
evaporation and convection, for our Titan model with latitudinal heat transport (a, solid 
line) exceed radiative fluxes in a model without heat transport and forced by the global- 
mean insolation (thin gray line) by a factor of six. In addition, there is a compensation 
between dry static energy fluxes into the surface (dotted line) and evaporative cooling 
(dashed line) which further enhances the methane cycle. In contrast, Earth's surface 
energetics (b) is dominated by direct solar heating, as is indicated by the net surface 
radiative imbalance in a model calculation without horizontal heat transport (thin gray 
line) as compared to the imbalance in NCEP reanalysis (solid line). Latitudinal transport 
accounts for a ~30% enhancement in the equatorial radiative imbalance. Both sensible 
(dotted line) and evaporative (dashed line) fluxes act to cool Earth's surface, except in 
very high latitudes. 
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Fig. 4. — The zonal-mean surface energetics of our model is dominated by a local balance 
between evaporative cooling (left) and radiative heating (right), with a relatively small role 
for sensible heat fluxes (center). Local, instantaneous surface radiative imbalance peaks at 
the summer pole, exactly when the top-of-atmosphere imbalance peaks (Figure [T]). 
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